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ABSTRACT 

Local  and  average  Nusselt  numbers  have  been  experimentally 
determined  around  the  circumference  of  a  uniformly  heated 
right  circular  cylinder  placed  near  a  plane  surface  in  a 
cross  flow  of  air.   Experiments  were  carried  out  for  Reynolds 
num'bers  ranging  from  16,000  to  153  >  000,  three  plate  lengths 
and  varying  cylinder  to  plate  spacings.   The  variations  of  the 
average  Nusselt  number  as  a  function  of  gap  spacing  was  found 
to  exhibit  an  interesting  and  unexpected  trend.   The  average 
Nusselt  number  exhibited  a  minimum  value  when  the  cylinder 
was  resting  on  the  plate  and  attained  a  maximum  value  at  a  gap 
spacing  of  approximately  one  to  two  cylinder  radii.   Flow 
visualization  studies  together  with  measurements  of  the  vortex 
shedding  frequency  in  the  cylinder  wake  indicate  that  when  the 
plate  is  within  approximately  one  cylinder  radius  of  the 
cylinder,  it  inhibits  the  formation  of  the  regular  vortex 
pattern  that  forms  in  the  wake  of  a  free  cylinder.   This  has 
a  strong  influence  on  the  base  Nusselt  number  which,  in  turn, 
exerts  a  controlling  influence  on  the  magnitude  of  the  average 
Nusselt  number. 
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I.   INTRODUCTION 

An  experimental  investigation  was  carried  out  to  examine 
the  influence  of  wall  proximity  effects  on  the  heat  transfer 
characteristics  of  a  uniformly  heated  cylinder  placed  in  a  cross- 
flow  of  air.   Recently,  McComas  /~1_7  initiated  an  investigation 
of  the  heat  transfer  characteristics  of  a  uniformly  heated 
cylinder  attached  to  a  plate  and  Gnerlich  /~2_7  continued  the 
same  study  varying  the  cylinder  to  plate  spacings.   The  present 
work  is  an  extensive  follow-on  study  of  the  work  of  McComas 
and  Gnerlich  differing  in  that  a  range  of  Reynolds  numbers  has 
been  investigated  as  well  as  three  different  plate  lengths. 
The  present  experimental  work  was  conducted  in  the  Mechanical 
Engineering  Department  wind  tunnel  (described  in  detail  in 
Section  III)  in  Building  500  at  the  Naval  Postgraduate  School. 

Local  Nusselt  numbers  v/ere  obtained  around  the  circumference 
of  the  cylinder  at  Reynolds  numbers  of  16,000,  90,000  and 
153,000  without  a  plate  as  well  as  with  a  plate.   In  the 
case  with  a  plate,  the  cylinder  to  plate  spacings  were  varied 
from  d/r  =  0  to  d/r  =5-33  with  L/D  =  k   where  d  =  gap  width, 
r  =  cylinder  radius,  D  =  cylinder  diameter  and  L  =  length  of 
plate  from  the  centerline  of  the  cylinder  to  the  leading  edge 
of  the  plate  (Figure  l) .   Experiments  were  also  conducted  at 
a  Reynolds  number  of  153,000  with  L/D  =  2  and  8  in  order  to 
see  the  effect  of  the  length  of  the  plate's  leading  edge  on 
the  heat  transfer  results. 
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Before  proceeding  further,  a  brief  qualitative  description 
of  flow  past  a  cylinder  will  help  to  explain  and  interpret  the 
results  obtained  from  this  experimental  study.   Since  the 
present  experiments  examined  Reynolds  numbers  ranging  from 
16,000  to  153,000,  only  the  subcritical  and  critical  flow 
regimes  will  be  described  here.   See  Achenbach  /~3_7  for  a 
complete  discussion  of  the  various  flow  regimes. 

In  subcritical  flow,  a  laminar  hydrodynamic  boundary 
layer  grows  from  a  minimum  thickness  at  the  forward  stagnation 
point  on  the  cylinder  to  a  maximum  thickness  at  an  angular 
location  of  approximately  80-85  degrees.   At  this  location, 
the  kinetic  energy  of  the  moving  fluid  is  no  longer  sufficient 
to  overcome  the  adverse  pressure  gradient  present  on  the 
cylinder's  surface  and  the  laminar  boundary  layer  separates. 
Upon  separation  a  turbulent  wake  then  develops  on  the  rear 
side  of  the  cylinder. 

In  critical  flow,  the  laminar  boundary  layer  again  grows 
from  the  forward  stagnation  point  to  the  laminar  separation 
point  which  occurs  at  approximately  80-85  degrees.   In  contrast 
to  the  subcritical  flow  case,  a  transition  occurs  from  laminar 
to  turbulent  flow.   A  turbulent  boundary  layer  then  reattaches 
itself  to  the  cylinder.   This  reattached  turbulent  layer 
eventually  separates  at  a  position  downstream  of  the  laminar 
separation  point.   The  region  between  the  laminar  separation 
point  and  the  point  of  reattachment  of  the  turbulent  boundary 
layer  is  known  as  a  "separation  bubble." 
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The  results  of  the  experiments  with  the  cylinder  attached 
to  the  flat  plate  indicate  a  critical  flow  pattern  exists 
on  the  top  half  of  the  cylinder  for  all  Reynolds  numbers 
investigated.   This  flow  pattern  changed  from  critical  to 
sub-critical  over  the  top  as  the  cylinder  to  plate  spacing 
was  increased  to  approximately  d/r  =  0.1  to  0.2.   The  flow 
over  the  bottom  half  of  the  cylinder  revealed  critical  flow 
immediately  after  a  space  existed.   Flow  over  the  bottom 
turned  from  critical  to  subcritical  in  the  interval  of 
0.8  ^  d/r  ^1.2. 

Plots  of  the  average  Nusselt  number  versus  gap  spacing 
have  shown  that  the  minimum  average  Nusselt  number  occurs  when 
the  plate  is  attached  to  the  cylinder  and  the  maximum  average 
Nusselt  number  occurs  at  d/r  -  2.0. 

Qualitative  studies  /~^_7>    using  a  flow  visualization 
technique  in  the  water  tunnel  together  with  measurements  of 
the  vortex  shedding  frequency  in  the  wind  tunnel,  have  contri- 
buted a  plausible  explanation  about  the  trends  of  average 
Nusselt  number  curves. 

It  is  theorized  that  trends  in  the  average  Nusselt 
number  curves  are  a  direct  consequence  of  the  strength  and 
frequency  of  the  vortex  pattern  which  forms  in  the  wake  of 
the  cylinder. 
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II.   BACKGROUND 

The  present  investigation  of  the  heat  transfer  character- 
istics of  a  uniformly  heated  cylinder  placed  in  close  proximity 
to  a  plane  surface  in  an  air  stream  was  conducted  concurrently 
with  an  investigation  by  G'dktun  £^J   of  the  flow  character- 
istics of  a  cylinder  placed  near  a  plane  surface.   Information 
obtained  from  the  survey  of  the  heat  transfer  experiments 
proved  useful  in  explaining  the  nature  of  the  flow  field 
observed  by  Goktun.   In  the  same  manner,  Goktun' s  observations 
made  significant  contributions  to  a  better  understanding  of 
the  heat  transfer  characteristics  of  the  cylinder. 

Although  extensive  investigations  have  been  directed  at 
obtaining  a  better  understanding  of  the  heat  transfer  and  flow 
characteristics  of  "free"  cylinders  (for  extensive  references 
on  this  subject  see  Cooper,  Field  and  Meyer  £ ~  5£7)  »  data 
relating  to  the  heat  transfer  and  flow  characteristics  of 
the  cylinder-plate  geometry  under  investigation  were  not 
available  in  the  literature.   Such  information  would  be  of 
obvious  value  in  the  design  of  heated  piping  systems  placed 
near  walls  and  in  the  design  of  heat  exchangers. 

As  previously  mentioned,  the  present  investigation  is 
a  follow-on  study  of  an  investigation  initiated  by  McComas 
£l_7   scad    continued  by  Gnerlich  £ ~2_7 '.   Although  McComas 
experimentally  determined  local  and  average  heat  transfer 
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coefficients  around  a  uniformly  heated  cylinder  attached  to 
a  wall  and  Gnerlich  did  likewise  for  various  cylinder  to 
plate  spacings,  both  of  these  investigators  found  no  references 
to  theoretical  work  on  this  subject. 

Convective  heat  transfer  is  significantly  affected  by 
the  fluid  flow  characteristics  occurring  adjacent  to  a  solid 
boundary.   Since  forced-convection  heat  transfer  is  a  result 
of  the  movement  of  fluids  and  the  mixing  of  fluid  particles, 
the  mechanism  of  fluid  flow  must  be  understood  in  order  to 
understand  the  mechanism  of  heat  transfer  and  to  explain  the 
phenomena  occurring  during  the  heat  transfer  process.   Because 
of  this,  this  project  was  conducted  concurrently  with  Goktun's 
£~^_J   study  of  the  flow  field  near  the  cylinder-plate  geometry. 
In  explaining  the  heat  transfer  results  obtained  in  this 
work,  frequent  reference  will  be  made  to  Goktun's  concurrent 
flow  investigation. 
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III.   EXPERIMENTAL  APPARATUS 

1.   Wind  Tunnel 

An  open  circuit,  subsonic  wind  tunnel  was  used  for  the 
experiments.   The  tunnel  was  specially  designed  by  Kenney 
Engineering  Corp.,  of  Monrovia,  California  for  research 
quality  work .   The  entire  exterior  dimensions  are  given  in 
Figure  2.   The  tunnel  has  a  test  section  which  is  20  inches 
wide,  28  inches  high  at  the  entrance  and  expands  to  31  inches 
at  the  exit  to  provide  for  boundary  layer  growth.   It  is 
powered  by  a  75  HP,  1750  RPM  variable  speed  controlled  D.C. 
motor  which  provides  a  belt  drive  to  an  adjustable  pitch 
vaneaxial  fan.   The  fan  is  45  inches  in  diameter  and  is 
capable  of  supplying  air  at  10,000  to  70,000  cfm,  producing 
wind  velocities  ranging  from  10  mph  (1^.7  ft/sec)  to  175  rnph 
(256.7  ft/sec).   The  tunnel  is  equipped  with  six  graded  mesh 
screens  to  control  turbulence.   A  2k-   inch  water  manometer 
connected  in  to  the  stilling  chamber  and  test  section  may  be 
used  to  measure  velocity.   Utilizing  a  hot  wire  anemometer, 
a  turbulence  intensity  level  of  0.2$  was  found  in  the  wind 
tunnel  at  a  free  stream  velocity  of  98. 1  ft/sec.   A  mercury 
thermometer  was  installed  in  the  stilling  chamber  to  measure 
the  free  stream  temperature. 
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2.   Liquid  Crystals 

Microencapsulated  cholesteric  liquid  crystals  were  used  to 
obtain  both  qualitative  and  quantitative  heat  transfer  and 
fluid  flow  information  on  the  uniformly  heated  cylinder.   The 
liquid  crystals  provide  a  visual  display  of  temperature  over 
discrete,  reproducible  temperature  ranges  (called  the  "event" 
temperature  range)  by  means  of  a  change  of  color  (see  Cooper, 
Field  and  Meyer  /~5J7) •      The  visible  color  spectrum,  which 
falls  within  a  certain  temperature  band,  has  an  accuracy  of 
0.2  F.   All  experiments  in  this  investigation  were  conducted 
using  a  liquid  crystal  designated  as  NCR  S-43  having  an  event 
temperature  range  of 

Red   :   109.8°F   (^-3.20C) 

Green:   111.0°F   (^3»9°C) 

Blue  s   112. 3°F   (*J4.6°C) 

3-   Ternsheet 

An  electrically  resistive  carbon  impregnated  paper  known 
as  "Ternsheet"  was  used  as  the  heating  element  in  these  studies. 
Ternsheet  is  a  thin,  highly  flexible,  fibrous  material  contain- 
ing no  wires  or  ribbons  and  therefore  is  free  of  localized  hot 
spots,  burnout,  and  breakage  problems.   The  nominal  thickness 
is  0.1  cm  and  the  electrical  resistivity  is  about  2.5  ohm-cm. 
Using  the  Ternsheet  it  is  possible  to  heat  large  areas  within 
less  than  2fo   deviation  from  point  to  point. 
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k.      Cylinder 

The  test  cylinder  utilized  in  obtaining  the  heat  transfer 
data  was  the  one  used  by  Gnerlich  /~2j7.   It  was  made  of  three 
cylinders  attached  end  to  end  (see  Figure  3)  having  a  total 
length  of  20  inches,  corresponding  to  the  test  section  width, 
and  a  diameter  of  3  inches.   The  two  end  sections  were  hollow 
wood  cylinders  connected  by  a  wooden  brace,  curved  to  match 
the  cylinder  contour.   The  central  cylinder,  henceforth 
designated  as  the  test  section,  was  formed  from  a  sheet  of 
Temsheet.   It  was  tightly  packed  with  glass  wool  which  served 
the  dual  purpose  of  providing  insulation  against  possible 
internal  heat  losses  due  to  conduction  arid  free  convection 
and  also  provided  the  paper  cylinder  with  added  resistance 
to  deformation.   The  circumference  of  the  Temsheet  cylinder 
was  marked  at  five  degree  angular  intervals  in  two  locations 
2  inches  apart  for  accurate  location  of  the  liquid  crystal 
isotherms  during  the  experiment.   Aluminum  foil  conducting 
tape  served  as  electrodes.   The  aluminum  tape  was  circumfer- 
entially  attached  on  the  inner  surface  of  the  Temsheet.   Silver 
conducting  paint  was  applied  to  the  interface  between  the  tape 
and  Temsheet  to  insure  good  electrical  contact  was  maintained 
between  the  two  surfaces.   Two  lead  wires  were  connected  to 
the  aluminum  foil  through  one  end  of  the  cylinder.   The  lead 
wires  were  attached  to  a  regulated  DC  power  supply  (model: 
LAMBDA  LK  3^5A  FM) .   The  resistance  of  the  Temsheet  was 
measured  by  a  digital  multimeter  (model:   FLUKE  8100A) . 
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5.   Plates 

The  plates  used  in  the  experiments  were  constructed  of 
plexiglass  and  were  0.5  inch  thick,  20  inches  wide  and  had 
total  lengths  of  ^0 ,  46  and  58  inches,  respectively.   The 
plates  were  installed  horizontally  in  the  wind  tunnel  in  such 
a  way  that  the  leading  edge  of  the  plates  extended  2,  4  and 
8  cylinder  diameters  (i.e.  6,  12  and  2^  inches,  respectively) 
ahead  of  the  cylinder  centerline .   The  plates  were  beveled 
on  the  under  side  of  the  leading  edge  with  an  angle  of  12 
degrees.   The  plates  were  capable  of  being  located  at  any 
position  relative  to  the  cylinder  (Figure  k),   from  attached  to 
a  maximum  distance  of  8  inches  away  (i.e.  d/r  =  5*33) •   The 
plates  were  moved  vertically  by  means  of  a  threaded  rod 
extending  through  the  tunnel  floor,  and  a  rear  brace  which 
locked  in  position  by  means  of  a  split  nut. 


23 


.a 
E 
cd 
</> 
(/> 

(T3 

S- 
CD 

■o 

>> 

CJ 

-a 

c: 
ea 

CD 
fO 
Q. 
CD 

-t-> 


Q. 
fO 

s- 
a> 
o 
-t-> 
o 


CD 

i- 

CD 


2k 


IV.   EXPERIMENTAL  PROCEDURE 

The  present  investigation  was  conducted  using  the  same 
test  facilities  as  employed  by  Gnerlich  /~2_7.   Gnerlich, 
however,  performed  experiments  only  for  a  Reynolds  number  of 
153.000  (based  on  the  cylinder  diameter)  and  a  plate  leading 
edge  of  12  inches.   In  the  present  work,  experiments  were 
conducted  at  Reynolds  numbers  of  16,000,  90,000  as  well  as 
153.000  and  the  plate  leading  edge  was  varied  from  6  to  2^ 
inches . 

To  initiate  the  experiment,  a  sufficient  voltage  (from 
12  to  15  volts)  was  applied  to  the  cylinder  to  insure  that  the 
liquid  crystals  reached  their  event  temperature  and  changed 
color  over  the  entire  surface  of  the  cylinder.   This  was  done 
to  check  for  hot  and  cold  spots  which  would  indicate  improper 
packing  of  the  glass  wool  inside  the  cylinder.   The  wind 
tunnel  was  then  started  and  set  to  a  predetermined  velocity 
at  which  it  was  maintained  throughout  the  run.   The  applied 
voltage  was  adjusted  until  the  first  faint  red  color  trace 
was  observed.   The  power  supply  was  then  turned  off  and  the 
multimeter  was  switched  from  DC  voltmeter  to  ohmmeter.   The 
same  procedure  was  also  followed  at  the  end  of  the  experiment. 
In  this  manner,  an  accurate  value  of  the  paper  resistance, 
which  was  taken  to  be  the  arithmetic  mean  of  initial  and 
final  resistance  values,  was  determined.   The  multimeter  was 
then  switched  back  to  the  DC  voltmeter  and  power  was  re-applied 
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to  the  test  section.   A  steady  state  condition  was  usually 
attained  in  from  five  to  ten  minutes  depending  upon  the 
wind  velocity.   Generally,  at  least  ^0  minutes  were  spent 
to  take  the  first  data  point  in  order  to  "be  extremely  certain 
about  steady  state  conditions. 

The  first  color  band  observed  on  the  cylinder  was  at  the 
location  of  the  value  of  the  lowest  local  Nusselt  number. 
The  voltage  was  then  increased  such  that  the  liquid  crystal 
color  bands  would  appear  at  various  angular  positions  around 
the  circumference  of  the  cylinder.   The  temperature  obtained 
in  this  manner  could  be  compared  to  the  air  stream  temperature 
which  was  obtained  from  a  thermometer  mounted  in  the  stilling 
chamber.   Using  this  temperature  difference  and  power  dissi- 
pation, the  local  Nusselt  number  based  on  the  cylinder  diameter 
could  be  obtained  around  the  cylinder  (see  Appendix  A). 

Due  to  the  method  of  construction,  there  was  a  seam  along 
the  length  of  the  test  cylinder.   Measurements  within  30 
degrees  of  a  sector  of  the  seam  were  considered  unreliable. 
The  majority  of  Nusselt  number  data  were  collected  with  the 
seam  located  at  180  degrees.   The  Nusselt  number  data  in  the 
vicinity  of  .180  degrees  were  taken  by  turning  the  cylinder 
such  that  the  seam  was  facing  forward  at  0  degrees.   Several 
other  seam  positions  were  also  checked  but  it  was  ultimately 
decided  that  the  0  and  180  degrees  locations  were  the  best 
for  the  seam  position.   In  all  experiments,  the  location  of 
the  liquid  crystal  color  transition  points  were  easily 
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distinguishable  around  the  circumference  of  cylinder  except 
in  the  region  of  0  to  45  degrees  and  160  to  200  degrees. 
Especially  in  the  forward  stagnation  region  the  temperature 
gradient  was  very  shallow  and  the  whole  region  appeared 
isothermal.   In  the  vicinity  of  the  rear  stagnation  point, 
the  temperature  gradient  was  indistinct  due  to  the  turbulent 
behavior  of  the  flow  in  the  cylinder  wake . 
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V.   RESULTS 

The  results  shown  in  this  section  were  obtained  from 

the  heat  transfer  experiments  discussed  in  Section  VI •   The 

results  are  classified  in  two  groups: 

A.  LOCAL  NUSSELT  NUMBERS,  Nu 

a.  Typical  data  and  a  sample  calculation  for  the  heat 
transfer  coefficient  and  the  local  Nusselt  number  are  shown 
in  Appendix  A . 

b.  Figures  5  through  12  show  the  variation  of  the  local 
Nusselt  number  with  gap  size  at  a  Reynolds  number  of  153,000 
and  a  plate  leading  length  of  12  inches  (L/D  =  k  ) .      These 
plots  also  show  the  variation  of  the  pressure  coefficient 

as  a  function  of  gap  size.   All  pressure  data  v/ere  taken 
from  Goktun's  /~^J7   results  and  are  informative  in  that  they 
show  the  relationship  between  surface  pressure  and  heat 
transfer  characteristics. 

c.  Additional  plots  of  the  Nusselt  number  versus  gap 
size  for  Reynolds  numbers  of  16,000,  90,000  (with  L/D  =  l±   ) 
and  153,000  (with  L/D  =  2  and  8  )  can  be  found  in  Appendix  D, 
Figures  2^  through  S5>      Local  Nusselt  number  distributions  on 
the  surface  of  a  free  cylinder  are  also  contained  in  Appendix  D, 
Figures  56  through  58. 

B.  AVERAGE  NUSSELT  NUMBERS,  Nu 

a.   A  graph  of  Wu  versus  d/r  (the  cylinder  to  plate 
spacings)  for  Reynolds  numbers  of  16,000,  90,000  and  153 , 000 
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and  L/D  =  4.0  is  shown  in  Figure  13 .   A  comparison  is  made 
with  Gnerlich's  Z~2_7  results  on  the  same  graph. 

b.  A  graph  of  Nu  versus  d/r  for  Reynolds  number  of 
153,000  and  L/D  =  2,  4  and  8  is  shown  in  Figure  Ik. 

c.  A  graph  of  normalized  drag  coefficient  (C,/C,   ), 
Strouhal  number  (S/Sog  )  and  average  Nusselt  number  (Nu/Nu^,  ) 
versus  d/r  is  shown  in  Figure  15.   The  "co"  subscript  refers 
to  values  obtained  with  no  plate  in  the  tunnel. 

d.  A  graph  of  normalized  Strouhal  number  and  local 
base  Nusselt  number  (Nu,  se/Nuco  ^ase)  versus  d/r  is  shown 
in  Figure  16. 
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Figure  5-   Local  Nusselt  number  and  pressure  coefficient  on  the 
surface  of  a  cylinder  placed  near  a  plane  surface  for 
Re  =  153,000,  L/D  =  k,    d/r  =0.0. 
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Figure  6.   Local  Nusselt  number  and  pressure  coefficient  on 

the  surface  of  a  cylinder  placed  near  a  plane  surface 

for 

Re  =  15   v00,  L/D  =  k,    d/r  =?  0.25- 
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Figure  7.   Local  Nusselt  number  and  pressure  coefficient  on  the 
surface  of  a  cylinder  placed  near  a  plane  surface  for 
Re  =  153,000,  L/D  =  k,    d/r  =0.5. 
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Figure  8.   Local  Nusselt  number  and  pressure  coefficient  on  the 
surface  of  a  cylinder  placed  near  a  plane  surface  for 
Re  =  153,000,  L/D  =  kt    d/r  =1.0. 
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Figure  9«   Local  Nusselt  number  and  pressure  coefficient  on  the 
surface  of  a  cylinder  placed  near  a  plane  surface  for 
Re  =  153,000,  L/D  =  4,  d/r  -  2.0. 
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Figure  10.  Local  Nusselt  number  and  pressure  coefficient  on  the 
surface  of  a  cylinder  placed  near  a  plane  surface  for 
Re  =  153.000,  L/D  =  4,  d/r  =3-0. 
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Figure  11.   Local  Nusselt  num"ber  and  pressure  coefficient  on  the 
surface  of  a  cylinder  placed  near  a  plane  surface 
for  Re  =  153,000,  L/D  =  k,    d/r  =  *K  0  . 
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Figure  12.   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for   Re  =  153,000,  L/D  =  k,    d/r  -  5-33. 
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VI.   DISCUSSION 

Because  of  the  relationship  between  the  flow  characteristics 
(i.e.  pressure  distribution  around  the  cylinder,  pressure  drag 
and  vortex  dynamics)  and  the  heat  transfer  characteristics  of 
the  cylinder,  the  parallel  flow  study  of  Goktun  £~kj   will 
be  referred  to  in  some  parts  of  the  following  discussion. 
Goktun  obtained  his  results  using  an  unheated  aluminum  cylinder 
placed  in  proximity  to  the  plane  surface. 

The  general  nature  and  trends  of  the  results  obtained 
from  the  present  investigation  can  be  discussed  in  two 
groups: 

A.   TRENDS  IN  LOCAL  NUSSELT  NUMBERS,  Nu 

At  d/r  =  0 

The  initial  heat  transfer  measurements  were  made  with  the 
uniformly  heated  cylinder  attached  to  the  plate  at  a  Reynolds 
number  of  153 > 000  with  L/D  =  h .      The  data  obtained  in  this 
arrangement  agreed  within  the  estimated  experimental  uncertainty 
of  ±5fo   (See  Appendix  B)  with  the  results  of  McComas  /^\J 
and  Gnerlich  £~ 2j .      A  local  maximum  and  minumum  value  of 
the  Nusselt  number  were  once  again  observed  near  the  point 
of  attachment  of  the  plate  and  cylinder  and  are  indicative  of 
a  large  trapped  vortex  as  proposed  by  McComas  Z~l_7«   ^n  "^e 
flow  visualization  experiment,  at  d/r  =  0,  Goktun  /~^_7  also 
observed  that  a  series  of  small,  circular,  three  dimensional 
vortices  were  also  present  in  the  upstream  cavity  as  well  as 
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a  large  turbulent  vortex.   The  large  turbulent  vortex  occurs 
within  the  separated  region  where  the  flow  separates  from  the 
flat  plate  arid  reattaches  to  the  cylinder.   The  point  of  flow 
reattachment  occurs  at  approximately  1^0  to  170  degrees  from 
the  point  of  attachment  of  the  plate  and  cylinder.   The 
existence  of  this  large  turbulent  vortex  produces  the  observed 
local  maximum  and  minimum  values  in  the  Nusselt  number  curve. 

The  existence  of  a  separation  bubble  on  the  top  of  the 
cylinder  was  the  obvious  mark  of  critical  flow.   A  laminar 
boundary  layer  started  to  grow  on  the  forward  half  of  the 
cylinder  and  separated  between  80-85  degrees  and  a  subsequent 
transition  to  turbulent  flow  occurred  resulting  in  the 
reattachment  of  a  turbulent  boundary  layer  at  about  100  -105  • 
The  reattached  turbulent  layer  ultimately  separated  at 
approximately  130-135  degrees. 

The  decrease  in  the  local  Nusselt  number  (see  Figure  5) 
at  the  aft  end  of  the  cylinder  was  unexpected.   The  flow 
visualization  studies  clearly  showed,  however,  that  when  the 
cylinder  is  resting  on  the  plate  a  large  wake  with  very  little 
fluid  motion  develops  on  the  downstream  side  of  the  cylinder. 
The  absence  of  fluid  motion  therefore  explains  the  decrease 
in  Nusselt  number  on  the  back  of  the  cylinder. 

At  d/r  =  0.25  to  d/r  =5.33 

Further  experimental  results  at  a  Reynolds  number  of 
153,000  were  obtained  with  the  plate  as  far  removed  from 
the  cylinder  as  possible,  d/r  =  5-33' 
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Examination  of  the  polar  plots  in  Figure  6  to  12  for 
significant  trends  or  unusual  nature  show  two  items  of 
interest: 

1.  A  shift  in  the  location  of  the  forward  stagnation 
point,  rear  stagnation  point  and  the  laminar  separation  points 
( see  Figure  17) . 

Starting  from  d/r  =0.25  and  moving  the  plate  to 
d/r  =  2.0,  it  appears  as  though  an  imaginary  axis  initially 
passing  through  3^0  and  160  degrees  is  pivoted  about  the 
origin.   As  this  axis  turns  clockwise  through  20  degrees  of 
arc  to  360  degrees,  the  angular  position  of  the  stagnation 
and  separation  points  changed  correspondingly  up  to  d/r  =  2, 
and  then  reached  a  constant  behaviour  (see  Figure  6  to  9). 

This  observed  shifting  of  stagnation  points  was 
definitely  not  affected  by  the  variation  of  the  length  of 
the  plate's  leading  edge  and/or  the  Reynolds  number  as  shown 
in  Appendix  D,  in  Figures  ^0  to  5$. 

2.  The  transition  from  critical  to  subcritical  flow. 
The  transition  from  critical  to  subcritical  flow  was 

investigated  in  order  to  find  the  d/r  ratio  associated  with 
the  transition  on  the  top  and  bottom  of  the  cylinder.   At 
a  Reynold's  number  of  153 » 000  and  L/D  =  k,    the  flow  was 
critical  over  the  top  of  the  cylinder  with  the  plate  attached 
and  was  subcritical  over  the  top  at  d/r  =0.25.   By  alternately 
varying  the  plate  to  cylinder  spacing,  carefully  adjusting  the 
voltage  and  noting  the  liquid  crystal  display,  the  transition 
point  from  critical  to  subcritical  flow  was  precisely  found 
at  d/r  =0.1. 
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A  transition  from  critical  to  subcritical  flow  over 
the  bottom  of  the  cylinder  was  found  at  approximately 
d/r  =1.2.   In  this  case,  although  the  pattern  of  the 
separation  bubble  displayed  by  the  liquid  crystal  band  was 
less  pronounced  than  over  the  top,  it  was  clearly  observed 
as  far  as  d/r  =  1.2. 

As  also  noted  by  Gnerlich  /~2_7,  flow  over  both  the 
top  and  bottom  of  the  cylinder  are  definitely  critical  for 
the  stated  d/r  ranges  but  there  is  a  significant  difference 
in  the  degree  of  criticality  on  the  top  and  bottom  of  the 
cylinder  as  far  as  the  magnitudes  of  the  local  Nusselt  numbers 
are  concerned.   This  is  seen  in  Figures  5  through  8. 

Further  study  was  also  conducted  in  order  to  see  the 
influence  of  the  length  of  the  plate's  leading  edge  and 
Reynolds  numbers.   The  results  are  shown  in  Table  1,   Table  1 
shows  the  transition  points  from  critical  to  subcritical  flow 
that  occurred  over  the  top  and  bottom  of  the  cylinder  as  a 
function  of  Reynolds  number  and  L/D  ratios.   It  is  worthwhile 
to  note  all  points  fell  within  the  range  of  0.1  -  d/r  -0.2 
for  the  top.   The  transition  points  for  the  bottom  were 
found  in  the  interval  of  0.8  -  d/r  -  1.2.   Also  to  be  noted 
in  Table  1  is  that  the  transition  from  critical  to  subcritical 
flow  at  various  d/r  ratios  is  apparently  independent  of 
L/D  ratios  but  is  mildly  dependent  on  Reynolds  number.   As 
the  Reynolds  number  increased,  the  d/r  ratios  for  transition 
on  the  top  of  cylinder  decrease  slightly;  on  the  bottom  of 
the  cylinder  the  behaviour  was  reversed. 
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Re 

1/ 

D  -  2 

L/D  :  k 

L/D  :   8 

Top 
at  "d/r" 

Bottom 
at  "d/r" 

Top 
at  "d/r" 

Bottom 
at  "d/r" 

Top 
at  "d/r" 

Bottom 
at  "d/r" 

16000 

0.20 

1.0 

0.20 

0.9 

0.2 

0.8 

90000 

0.15 

0.90 

0.15 

0.9 

0.15 

0.8 

153000 

0.15 

1.2 

0.10 

1.2 

0.1 

1.2 

Table  1.   Transition  from  critical  to  suberitical  flow  as  a 
function  of  Reynolds  number,  d/r  and  L/D  ratios. 

B.   TRENDS  IN  AVERAGE  NUSSELT  NUMBER,  Nu 

In  order  to  gain  further  information  about  the  cylinder 
heat  transfer  characteristics,  the  average  Nusselt  number, 
Nu  was  computed  from  the  formula 


Nu  =  ^ 


2TT 
/     Nu-de 


The  integration  was  carried  out  on  an  IBM-360  computer  using 
Simpson's  Rule  of  Quadrature  with  a  step  size  of  5  degrees. 

Figure  13  shows  the  variation  of  the  average  Nusselt 
number  as  a  function  of  gap  spacing  for  three  Reynolds 
number  and  L/D  =  4 . 

A  comparison  of  the  results  obtained  at  Re  =  153 . 000  and 
L/D  =  k,    v/ith  those  obtained  by  Gnerlich  indicates  that  the 
trends  are  essentially  the  same.   The  decrease  in  Nu  at 
d/r  =  0.5  in  Gnerlich5 s  results  was  not  observed  in  the 
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present  experiments  but,  considering  the  magnitude  of  the 
decrease  and  the  Sf°   experimental  uncertainty  involved  in 
both  experiments,  it  is  not  felt  to  be  significant.   The 
maximum  Nu  was  observed  at  d/r  =  2.0  and  it  is  readily 
apparent  that  the  presence  of  the  plate  radically  altered 
the  magnitude  of  Nu  in  the  range  0  -  d/r  -  2.   Beyond  a 
d/r  of  2.0,  Nu  approached  a  nearly  constant  value. 

It  is  important  to  emphasize  that  the  trends  of  the  Nu 
versus  d/r  curves  are  nearly  the  same  for  all  Reynolds 
numbers . 

In  order  to  see  the  effect  of  the  length  of  the  plate's 
leading  edge,  results  were  obtained  for  a  Reynolds  number  of 
153,000  for  L/D  ratios  of  2,  k   and  8.   The  results,  which  are 
shown  in  Figure  14,  clearly  indicated  that  there  is  no 
significant  effect  on  the  trends  in  the  Nu  curve  due  to 
changing  the  L/D  ratios.   Goktun  /~^_7  also  noted  that 
changing  the  L/D  ratio  had  no  significant  influence  on  the 
drag  and  lift  characteristics  of  the  cylinder. 

For  the  sake  of  simplicity,  the  flow  characteristics 
(pressure  drag  and  vortex  shedding  frequency)  and  the  heat 
transfer  characteristics  (average  and  base  local  Nusselt 
number)  are  combined  and  plotted  in  normalized  form  in 
Figures  15  and  16. 

Goktun' s  qualitative  flow  visualization  study  /~^_7  in 
the  water  tunnel  together  with  the  vortex  shedding  frequency 
measurements  in  the  wind  tunnel  indicated  that  vortex  dynamics 
control  the  drag  on  the  cylinder.   Figure  15  clearly  indicates 
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the  interrelation  betv/een  the  drag  and  heat  Transfer .   It  is 
seen  that  as  the  drag  coefficient  increases,  the  average  heat 
transfer  coefficient  also  increases.   It  is  therefore  concluded 
that  the  vortex  dynamics  are  responsible  for  the  trends 
observed  in  both  the  drag  and  average  Nusselt  number  curves. 
In  the  range  0<d/r<(l,  Goktun  found  that  the  natural  or 
regular  vortex  pattern  formed  in  the  wake  of  a  free  cylinder 
is  inhibited  by  the  presence  of  the  plate.   The  plate  signifi- 
cantly reduces  the  vortex  shedding  frequency.   This  in  turn  has 
an  effect  on  the  base  pressure  (see  Figures  5  through  8)  and 
base  Nusselt  number  (Figure  16) .   The  more  stagnant  the  wake, 
the  higher  the  base  pressure  and  the  lower  the  base  Nusselt 
number.   Therefore,  the  stagnant  wake  results  in  poor  heat 
transfer.   This  trend,  in  turn,  influences  the  drag  and 
average  Nusselt  number  (Figure  15) . 

Goktun  /T^J7   observed  that  strong  and  violent  vortices 
were  formed  behind  the  cylinder  when  d/r  was  approximately 
1.0  to  2.0.   These  vortices  transfer  high-momentum  cool  fluid 
from  outside  the  boundary  layer  in  toward  the  rear  surface 
of  the  cylinder.   The  cool  fluid  in  the  vortex  which  impacts 
on  the  hot  cylinder  surface  has  a  large  capacity  to  remove 
energy.   The  "scrubbing"  motion  of  these  vortices  building 
up  behind  the  cylinder  and  breaking  away  to  be  carried  down- 
stream in  the  v/ake,  results  in  high  rates  of  heat  transfer 
and  this  results  in  a  large  base  Nusselt  number.   This  trend, 
in  turn,  influences  both  the  drag  coefficient  and  the  average 
Nusselt  number. 
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As  the  plate  is  moved,  further  from  the  cylinder,  the 
drag  and  Nusselt  number  approach  their  "free"  cylinder 
values . 
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VII.   SUMMARY 

This  experimental  investigation  was  carried  out  to  gain 
an  understanding  of  the  heat  transfer  behavior  of  a  uniformly 
heated  right  circular  cylinder  placed  near  a  plane  surface. 
A  concurrent  study  by  G ok tun  f^J   was  also  carried  out  to 
study  the  lift  and  drag  characteristics  of  a  cylinder  in  the 
same  configuration.   Local  Nusselt  numbers  v/ere  measured  around 
the  cylinder  circumference  for  various  plate  to  cylinder 
spacings,  plate  lengths  and  Reynolds  numbers.   The  average 
Nusselt  numbers  were  obtained  from  the  local  Nusselt  number 
distributions  by  numerical  integration. 

An  unusual  trend  which  was  observed  in  the  plots  of  both 
the  average  Nusselt  number  and  overall  drag  versus  cylinder  to 
plate  separation  curves  lead  to  a  flow  visualization  experiment 
and  hotwire  anemometer  measurements  of  the  vortex  shedding 
frequency  in  the  cylinder  wake.   Based  on  these  observations 
and  measurements,  it  was  concluded  that  the  structure  of 
vortex  pattern  that  develops  in  the  wake  of  the  cylinder 
exerts  a  controlling  influence  on  the  shape  of  the  average 
Nusselt  number  and  drag  curves.   In  the  range  of  a  plate  to 
cylinder  spacing  of  one  cylinder  radius,  it  was  found  that 
the  presence  of  the  plate  inhibits  the  formation  cf  the 
regular  vortex  pattern  that  forms  in  the  wake  of  a  free 
cylinder.   This  significantly  affected  both  the  cylinder 
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base  pressure  and  local  Nusselt  number  distribution  which, 
in  turn,  influenced  the  drag  and  average  Nusselt  number 
curves. 
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APPENDIX  A 

HEAT  TRANSFER  DATA  REDUCTION 

The  procedure  followed  in  reducing  the  heat  transfer 
data  is  shown  "below.   In  this  manner  the  local  Nusselt 
numbers  at  various  locations  were  obtained  for  the  case 
of  a  uniformly  heated  (constant  heat  flux)  cylinder. 

HEAT  FLUX 

The  heat  flux  produced  by  Joulean  heating  is  given  by; 

qw  "   A'R 

where:   E   =  Voltage  applied  to  the  test  section,   /~volts_7 
R   =  Total  resistance  of  Temsheet,  ^ohms_7 

A   =  Total  heated  area  of  Temsheet  (  =0.315)  /"ft  J 

v2 
q   =  Rate  of  total  heat  transfer,  v 

n-  fir 


V2       0  ,nQ    BTU 

Conversions  1  5  -  3-^-lJ  0" 

n-  f1T  hr-ft^ 


Hence,   4W  =  ^~    4     =  »•««  "£     £*? 
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HEAT  TRANSFER  COEFFICIENT 

The  total  heat  transfer  coefficient  is  determined  by 
dividing  the  heat  flux  by  the  measured  local  temperature 
difference  (T  -  T^  )  : 


ht  = 

*w 

^w 

T  -  T^ 

AT 

ht  = 

h  +  h 
c    r 

h  = 

0 

4w 

AT      nr 

hr  =  a'3f-2  *(T  +  T-  )  (T2  +  T-  } 


where : 


BTU 


h,  =  total  heat  transfer  coefficient, 
X  hr-ft  -  °F 

T  =  cylinder  surface  temperature  at  angular 

location  e  ,    as  determined  by  liquid  crystals,     F 

Top  ~  air  stream  temperature,  F 

"RTTT 

h  =  convective  heat  transfer  coefficient, 
c 

h  =  radiation  heat  transfer  coefficient, 


o-  =  Stefan-Boltzman  constant  =  0.171^  x  10 

hr   -    fl;      -    °R^ 


(p-io"   ra-diation   exchange   factor  between   cylinder 
and   surroundings  =      0.9 


hr  -  ft2  - 

■  °F 

BTU 

hr  -  ft2  - 

■  °F 

-8     BTU 

5^ 


NUSSELT  NUMBER 


h  •  D 
Nu  =  — 


kf 


where:   D  =  cylinder  diameter  =  i-  f t 

kf  =  thermal  conductivity  of  air  evaluated  at 


m    i   m 

x     x°°  BTU 


film  temperature  T„  =  

2       '  hr  -  ft  -  °F 


NUMERICAL  EXAMPLE: 

A  sample  calculation  of  the  local  Nusselt  number  is 
shown  below  using  the  data  (see  page  5?   )  obtained  at  a 
Reynolds  number  of  153,000,  d/r  =  4.0,  L/D  =  8  and  angular 
location  of  0-6^   degrees. 

V  =  31  volts 

R  -  13-3 

T  -  111°F  -  571  °R 

TccP  70°F  =  530°R 

£T  =  T  -  TTO  =  4l°F 

9T-2  =  °-9 

0    =    0.1714  X  10"8   ggg,   0  U 

hr  -  ft     -  R 

1  • 
D  =  3"  =  T, 

4  =  10.842  I-   =  10.8^2  f2il  =  783.39    BTU  0 
^w  R  13.3  hr  _  fi^ 

h  =  _t_  a   783.39   =  19(1     BTU 


t    AT       41  hr  -  ft2  -  °F 
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hr  =  *'K-2    (T  +  T~  )(T2  +  T<*  ) 

=  (0.1714  X  10"8)  (0.9)  (571  +  530)  (5712  +  5302) 
=   1  BTU 


hr  -  ft2  -  °F 


h  =  h.  -  h  =19. 1-1=  18.1 ^~ 

hr  -  ft   -   F 

m  +  m 

00      111  +  70         o 
Tf  =  =  ■   '  -^  (        =    90.5  F 


k„  =  0.0154  (evaluated  at  film  temperature 

1  hr  -  ft  -  °F  of  90.5  °F) 


Nu  = 


hc  'D    _   (18.1)  (  5J. 


kf  0.0154 


Nu  =  29^   at  6 =  63  degrees. 
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L/D=   8.< 

DATA 

uoo"  98.1     ft/sec 
«  2.2       in-H20 

RT  -  1J.3    Ohm. 

», 

3     ■ 

153 

,000 

DATE:  23/  Sept./ic?75 

TIMS 

(volts; 

Jry.Te, 

kid    To 

to  0) 

hi 

nr 

hc 

Ku 

e 

COMMENTS 

1350 
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APPENDIX  B 

UNCERTAINTY  ANALYSIS 

Since  there  are  uncertainties  in  the  measured  quantities, 

an  uncertainty  analysis  was  performed  to  ascertain  the  degree 

of  uncertainty  in  the  experimental  results  at  three  ranges  of 

Reynolds  numbers  of  153,000,  90,000  and  16,000.   The  method 

of  Kline  and  McClintock  /~6_7  was  used. 

Being  conservative,  the  following  uncertainties  were 

estimated  in  the  measured  variables: 

voltage,   V  *  0.1  volts  (20:1) 

resistance,  R  t  O.A1-  ohms  (20:1) 

temperature,   T  ±  1 . 0  °F  (20:1) 

area,     A  t    0.001  ft2  (20:1) 

diameter,      D  t    0.001   ft  (20:1) 

wall  manometer,    H  1    0.05    in   -   H2o  (20:1) 

radiation  exchange 

factor        ff?   t  0.05  (20:1) 

The  uncertainties  associated  with  the  properties  of  air 
were  considered  to  be  negligible. 

The  uncertainty  analysis  is  performed  for  the  data 
obtained  at  Reynolds  numbers  of  153.000,  90,000  and  16,000 
as  follows; 
1 .   Illustration  of  procedure  of  uncertainty  analysis  for 

Re  =  153.000  ! 
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Total  Heat  Transfer  Coefficient,  h,  : 

As  shown  previously  in  Appendix  A,  the  total  heat 
transfer  coefficient  was  calculated  from  the  equation: 


2  2 

.   =  V  V 

nt  "   R  •  A  •  (T-Tco  )   "   R  .  A  •  AT 


The  uncertainty  in  the  heat  transfer  coefficient  is  then 
expressed  as: 


X     _    ./    (Z%     \\    (<*r    f,      (   "A    Y    ,     (__a_AT 


h.  W\V/  \R/  \A      /  \aT 

For  the  run  at  a  Reynolds  number  of  153.000,  L/D  =  8, 
d/r  --  4  and  an  angular  location  of  63  degrees: 
V    =31-0.1  volts 
R    =13.3-0.4  ohms 

AT    =  41  ±  1.0  °F  (  T  =  111°F  =  571°R.  ^  =  70°F  =  530°R) 
A    =  0.314  t  0.001  ft2 
H    =  2.20  t  o.05  in  -  H2o 

BTU 


%+   =19.1 


t  hr  -  ft^  -  °F 


h_t_  =  0.04   or  u.    =  t  0.8   BTU 


ht  t  hr  -  ft2  -  °F 

This  result  does  not  take  into  consideration  the 
uncertainty  in  the  angular  location  of  the  liquid  crystal 
band.   This  uncertainty  is  estimated  to  be  -  3  degrees. 
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The  uncertainty  in  the  total  heat  transfer  coefficient 
remained  relatively  constant  around  the  circumference  of  the 
cylinder.   However,  the  uncertainty  in  the  angular  location 
of  the  liquid  crystal  hand  varied.   An  estimated  maximum 
uncertainty  of  -  5  degrees  occurred  in  the  forward  and  rear 
stagnation  regions  and  a  minimum  of  -  1  degree  occurred  in  the 
regions  near  the  separation  and  reattachment  points. 


Radiation  Heat  Transfer  Coefficient,  h   : 

I,, 

The  equation  used  to  calculate   h   was: 

2  ,  m2 


hr  =  a-  •  J^_2  .(T  +  T^  )(:r  +  T^  ) 


1 


defining  T_   -~-   (  T  +  T^  ) 


r  1-2    m 


Since  a  is  constant,  it  is  not  considered  in  the  uncertainty 
analysis.   The  uncertainty  in  h   is  then  expressed  as: 


h 
r 


% 


+ 


BTU 


h—  =  0-06    0r   »h   =  °'°6   ._    „2 


r 


h  =  1.00  t  0.06 


hr  -  ft"  -   F 

BTU 


_L.  «   V  \-J  -KS      w      v     *-r  — 


hr  -  ft^  -  °p 
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Local  Nusselt  number,  Nu: 

The  local  Nusselt  number  was  calculated  from  the 

equation: 

h   -D 
Nu  = 


kf 


Therefore , 


Wl\T 


Nu 


Nu 


Note  that;   h  =  h .  -  h 
c    t    r 


wh 
c 


h 
c 


Since  the  uncertainties  in  the  properties  of  air  are  considered 
to  be  negligible  and  the  uncertainty  in  the  cylinder  diameter 
is  negligibly  small; 


2  v  2 


wh,  \       /  wh 


Nu     /  /   t    \     +   /  __r 


Nu     1/      h   J  V  h„ 

c  /      \   c 


=   0.0^4    or  „,„„  -  13 


NU  v/.w-x-r       ~X   UJ^ 

Nu  =  29^  -  13 
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Reynolds  Number,  Re  ; 

The  equation  used  to  calculate  Reynolds  number  v/as : 
U   -D 


Re  = 


\y 


Therefore 


Re 
Re 


Since  the  uncertainty  in  the  kinematic  viscosity  and  in 
the  cylinder  diameter  are  negligible: 

OV,         U 

Re       °° 


The  quantity  Ucc  is  proportional  to  the  square  root  of 
the  wall  manometer  reading,  hence: 


0). 


Re      1    \ 


Re      2    H 


wRe 

-   0.011   or  wD   =  2000 

Re 


Re 


Re  =  153,000   -  2000 


2.   The  uncertainty  analysis  for  the  run  at  Re  =  90 . 000: 

For  the  data  obtained  at  a  Reynolds  number  of  90,000, 
L/D  =  k,    d/r  =  2.0  and  6 =  67  degrees: 
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V   =  28.0  -  0.1   volts 

R   =13.0-0.4   ohms 

AT  -  38.6  t  1  °f   (  T  =  111°F  =  571°R,  ^  =  72. 4°F  =  532. 4°R) 

A   =  0.314  -  0.001  ft2 

H   =  O.76  -  0.05   in  -  Ho 

Similarly,  the  uncertainties  calculated  in  following 
variables  were: 

Wh 

-r^-   =0.04    or   to,   =  t  O.?   ^- 

ht  ht  hr  -  ft2  -  °F 

h.  =  16.9  -  0.7  ~- 

t  2    o 

hr  -  f1T  -   F 

■  =      0.06   or  to.    =   0.06 


hr  hr  hr  _  ft2  _  op 


BTU 


h   =  1.00  t  0.06  

r  .     o,2   o„ 

hr  -  ft   -   F 


to 

=  0.044   or  ukt  -  11 


Nu     -•—   -x   ^Nu 


Nu  =  259   -  11 


to 

=   .032       o)Rq=  3000 


Re       J  Re 

Re   =  90,000  t  3000 
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3-   The  uncertainty  analysis  for  the  run  at  Re  =  16,000  : 

For  the  data  obtained  at  a  Reynolds  number  of  16,000 
L/D  =  4,  d/r  -  2.0  and  0  =  65   degrees: 


V  =21.1   -  0.1   volts 

R  =13.4  t  0.4   ohms 

AT  -    43     t    1   °F    (    T=    111°F  =    571°R,    Tco  =    68°F   =    528°R) 

A  =    0.314      -    0.001    ft2 

The  following  uncertainties  were  calculated: 


0.039   or  03,   =  t  0.3     BTU 


j.  h 


"t  t  hr  _  ft2  _  oF 

ht  =  8.3  ±  0.3   BTU 

hr  -  ft2  -  °F 


wh 


,—  =  0.06    or  w,   =  t  0.06   ^TU 

r  r  hr  -  ft  -  °F 

h   =  1.00  t  0.06  mi 

hr  -  ft^  -  °F 


0JNu 
Nu 


=  0.04    or  ooNu  =  -  5 
Nu  =120   -  5 


For  this  particular  Reynolds  number  it  was  impossible  to 
measure  the  free  stream  velocity  using  the  wall  manometer. 
The  hot  wire  anemometer  was  therefore  used  to  measure  the 
velocity.   Thus,  the  uncertainty  estimated  in  U^  was: 
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U^  ±  2  ft/sec  (20:1) 


CO 


Re 


Re   "=  U^~-  =  -10-  =  °'Z      or  WRe  =   3°00 

Re  =  16,000   -  3000 

The  high  uncertainty  in  Uoo  and  consequently  in  the  Reynolds 
number  was  primarily  due  to  the  unstable  operation  of  the 
fan  motor  at  its  lowest  operating  point. 
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APPENDIX  C 
WIND  TUNNEL  VELOCITY  PROFILES 

For  the  three  ranges  of  Reynolds  numbers  of  16,000, 
90,000  and  153,000  and  with  L/D  =  k,    velocity  profiles  were 
obtained  in  the  tunnel  test  section  (see  Figure  18  through 
23).   Data  were  collected  at  d/r  =  0  and  d/r  =  5.33  (see  next 
page  for  data) . 

The  hot  wire  anemometer  was  used  to  probe  the  test 
section  at  various  positive  and  negative  Y  -  locations  at 
90  degrees  and  270  degrees,  measured  from  the  forward 
stagnation  point  on  the  cylinder.   The  velocity  profile  at 
the  test  section  entrance  v/as  also  measured  at  various 
locations  and  it  v/as  found  to  be  fairly  uniform. 

The  turbulence  intensity  in  the  test  section  was  obtained 
for  a  Reynolds  number  of  153*000  and  was  found  to  be  approxi- 
mately 0.2%. 
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LATA  for   the  velocity  profiles 
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Figure   18. 


Velocity  profile   at 
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Figure   19 . 


Velocity  profiles    at   the    top   and   bottom   of  the 
cylinder   for   d/r  =    5.33. 
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Figure  21. 


Velocity  profiles  at  the  top  and  bottom  of  the 
cylinder  for  d/r  =  5»33« 
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Figure  22.   Velocity  profile  at  the  top  of  the  cylinder  for 
d/r  =0.0. 
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Figure  23.   Velocity  profiles  at  the  top  and  bottom  of  the 
cylinder  for  d/r  =  5-33. 
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APPENDIX  D 

Additional  plots  of  the  Nusselt  number  versus  gap  size 
for  Reynolds  numbers  of  16,000,  90,000  and  153,000. 
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Nu 


Figure  2k.      Local  Nusselt  number  on  the  surface  of  a  cylinder 
placed  near  a  plane  surface  for 
Re  =  16,000,  L/D  =  k,      d/r  =0.0. 
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Nu 


Figure  25.   Local  Nusselt  number  on  the  surface  of  a  cylinder 
placed  near  a  plane  surface  for 
Re  =  16,000,   L/D  =  k,      d/r  =  0.25- 
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Figure  26.   Local  Nusselt  number  on  the  surface  of  a  cylinder 
placed  near  a  plane  surface  for 
Re  =  16,000,   L/D  =  k,      d/r  =0.5- 
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Nu 


Figure  27.   Local  Nusselt  number  on  the  surface  of  a  cylinder 
placed  near  a  plane  surface  for 
Re  =  16,000,   L/D  =  k,      d/r  =1.0. 


78 


Nu 


_U_UJ_L1_LJ 


Figure  28.   Local  Nusselt  number  on  the  surface  of  a  cylinder 
placed  near  a  plane  surface  for 
Re  =  16,000,   L/D  =  4,   d/r  =  2.0. 
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Figure  29.   Local  Nusselt  numoer  on  the  surface  of  a  cylinder 
placed  near  a  plane  surface  for 
Re  =  16,000,   L/D  =  k,      d/r  =3.0. 
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Figure  30 .   Local  Nusselt  number  on  the  surface  of  a  cylinder 
placed  near  a  plane  surface  for 
Re  =  16,000,   L/D  =  k,      d/r  =  k.Q, 
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Figure  31  •   Local  Nusselt  num"ber  on  the  surface  of  a  cylinder 
placed,  near  a  plane  surface  for 
Re  =  16,000,  L/D  =  kt      d/r  =  5-33. 
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Figure  32.   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  90,000,  L/D  =  k,    d/r  =0.0. 


83 


Nu   • 
Cp  x 


760 


Figure  33-   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  90,000,  L/D  =  k,    d/r  =  0.25. 
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Figure  Jk .      Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  90,000,  L/D  =  4,  d/r  =0.5- 
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Figure  35 •   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  90,000,  L/D  =  k,    d/r  =1.0. 
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Figure  36.   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for   -Re  =  90,000,  L/D  =  k,    d/r  =  2.0. 
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Figure  37  •   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  90,000,  L/D  =  k,    d/r  =  3-0. 
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Figure  38.   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  -  90,000,   L/D  =  4,  d/r  =  4.0. 
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Figure  39 «   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  90,000,  L/D  =  k,      d/r  =  5-33. 
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Figure  ko  .      Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  =  2,  d/r  =0.0. 
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Figure  kl.      Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  =  2,  d/r  =  0.25. 


92 


Nu   • 
Cp  x 


Oi.Z 


"260 


Figure  k2.      Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  =  2,  d/r  =0.5- 
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Figure  ^3 •   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  =  2,  d/r  =1.0. 
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Figure  kk .      Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153.000,  L/D  =  2,  d/r  =  2,). 
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Figure  45 .   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  =  2,  d/r  =3.0. 
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Figure  ^-6.   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  -  153.000,   L/D  =  2,  d/r  =  *K0. 
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Figure  k7 .      Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  =  2,  d/r  =  5-33. 
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Figure  ^-8.   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153.000,  L/D  =  8,  d/r  =0.0. 
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Figure  ^9-   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  =  8,  d/r  =  0.25. 
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Figure  50.   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  =  8,  d/r  =  0.5- 
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Figure  51  •   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  =  8,  d/r  =1.0. 
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Figure  52.   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  =  8,   d/r  =  2.0. 
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Figure  53  •   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,   L/D  =  8,  d/r  =3-0. 
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Figure  5^«   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  =  153,000,  L/D  -  8,  d/r  =4.0. 
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Figure  55  >      Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  cylinder  placed  near  a  plane 
surface  for  Re  -  153,000,  L/D  =  8,  d/r  =  5-33. 
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Figure  $6,      Local  Nusselt  number  on  the  surface  of  a  free 
cylinder  for  Re  =  16,000. 
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Figure  $7 .   Local  Nusselt  number  and  pressure  coefficient  on 
the  surface  of  a  free  cylinder  for  Re  =  90,000. 
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Figure  58.   Local  Nusselt  number  and  pressure  coefficient 

on  the  surface  of  a  free  cylinder  for  Re  =  153.000 
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